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Abstract 



We use lattice methods to evaluate from first principles the spectrum of hybrid mesons 
produced by gluonic excitations in quenched QCD with quark masses near the strange quark 
mass. For the spin-exotic mesons with J PC = 1 h , 0^ , and 2 H which are not present in 
+^ ■ the quark model, we determine the lightest state to be 1 ^ with mass of 2.0(2) GeV. 

One of the goals of quantitative studies of QCD is to determine the masses and properties of 
states which are not allowed in the simple quark model because they contain gluonic excitations. 
The prototype for such a state is the glueball and accurate lattice studies |2| have been able 
to pinpoint the relevant mass range for experimental study of the scalar glueball as around 1.6 
GeV. Another important area is the study of hybrid mesons which are qq mesons with gluonic 
excitation. Because the gluonic excitation can introduce angular momentum, the most clear cut 
signal for a hybrid meson is to search for J PC quantum numbers not allowed in the quark model. 
These include J PC = 1"+, 0+~, and 2+~, A gain lattice QCD is capable, from first principles, 
of establishing the masses of these states. 

In full QCD, the quark pair creation and anihilation processes will lead to the mixing of 
glueballs with qq mesons and it is thus appropriate to establish the glueball spectrum with 
these processes turned off - the quenched approximation. This then acts as a guide for mixing 
studies. For hybrid mesons, however, the quenched approximation already allows mixing between 
qq mesons and hybrid mesons if their J PC quantum numbers are non-exotic. This will be 
quite difficult to untangle in lattice studies, so we focus instead on the exotic J PC cases. Any 
mixing for these exotic J PC states would have to be with qqqq mesons and such mixing is 
turned off in the quenched approximation. Thus, in the quenched approximation, we are able 
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to determine unambiguously the exotic hybrid energy levels given by the full non-perturbative 
gluonic interaction. 

The spectrum of hybrid mesons has been explored in lattice QCD for many years [||, Q in 
the context of heavy quarks. In the quenched approximation, the lightest hybrid level is found 
to be degenerate and to contain spin-exotic mesons with J PC = 1 '", H , and 2 H . Since 
the quenched approximation does not reproduce exactly the experimental cc and bb spectra, the 
prescription of increasing the short distance component in the quenched approximation potential 
yields the lattice prediction || for the lightest hybrid meson excitation to be at 4.19(15) GeV for 
cc and 10.81(25) GeV for bb. These energy values lie above the open DD and BB thresholds. An 
alternative procedure, within the quenched approximation, is to focus || on the energy difference 
between the hybrid meson and the BB threshold. This suggests that the lowest hybrid level 
may lie below the threshold. If the hybrid states were to lie below these thresholds, then they 
would be narrow resonances and so much clearer to detect experimentally. These lattice results 
are for the quenched approximation and do not include mixing with 4-quark states. They also 
they do not include spin-orbit effects which may be significant for quarks of finite mass. This is 
one motivation for exploring hybrid levels on a lattice with realistic quark masses. 

There have been many experimental claims of signals for J PC exotic mesons in light quark 
spectroscopy (by light we mean composed of u, d, s quarks). None of these claims are over- 
whelming H 0] and guidance from QCD is sorely needed, as it was for glueball exploration. 
With present computational resources, it is possible to explore this and here we present the 
results of the first systematic study. Details of our lattice techniques and our preliminary results 
from much lower statistics have been published earlier S. Another lattice collaboration is also 
engaged in a study of this area || . 

Spin-exotic hybrid mesons have gluonic components with non-trivial angular momentum. 
This can be modelled on a lattice by joining the quark and antiquark by a colour flux path 
(product of links) which is not straight. Guidance in designing suitable paths comes from the 
heavy quark studies where a U-shaped path (actually a difference of the form n — U ) between 
static quarks was found ||, ||| to yield the lightest hybrid excitation. We construct hybrid meson 
operators for light quarks by joining the light quarks by colour flux in that U-shaped orientation. 
This involves an extra quark propagator inversion, but we are able to study a range of different 
hybrid mesons using this extra inversion. Effectively the technique is to employ a 'white' source 
on the lattice and to separate the hybrid state of interest by choosing an appropriate operator 
for the sink ||. 

Details of the lattice evaluation and operator construction are given in ref. || and we sum- 
marise them here. We used a 16 3 x 48 lattice at = 6.0 with the Wilson gauge action for SU(3) 
colour fields. This corresponds to a lattice spacing of approximately 0.1 fm. The quark prop- 
agators were evaluated at hopping parameter K = 0.137 using an SW-clover fermionic action 
with coefficient c = 1.4785. For our U-shaped source we took a separation of 6 lattice units and 
we evaluated 350 local propagators from sources at (0, 0, 0, 0) and at (0, 0, 6, 0). For the sink at 
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meson 



mass 



M/My M s /M((f>) 



multiplet 



Ma 



P 
a 

&2 



1-+, 3-+ 
0+-, 4+- 
2+" 4+- 
2+-, 3+- 



0.95(7) 
1.05(7) 
1.26(13) 
1.12(10) 



1.76(13) 
1.94(13) 
2.33(24) 
2.07(18) 



1.95(13) 
2.16(13) 
2.62(24) 
2.32(18) 



Table 1: The masses of the hybrid mesons. The J PC values are the lowest two values allowed 
by the lattice cubic symmetry. The masses are given in lattice units and as a ratio to the vector 
meson mass my. The mass ratio appropriate to ss states is given in the last column. 

time t, we used U-shapes of size 1, 3 and 6 lattice units. 

Using the same operator at source and sink (namely with size 6) gives an effective mass 
which is an upper bound on the ground state mass. In order to estimate this ground state mass, 
we evaluate the correlation between our source operator and all three sink operators since this 
data set allows us to estimate the excited state contribution and so extract the ground state 
mass by making stable two state fits to these correlations. The effective mass values for p are 
shown in the Figure. Given that the signal to noise ratio deteriorates rapidly with increasing 
time separation t, it is vital to use the small t data to extract the ground state mass. We used 
correlated fits to as wide a range of t as possible consistent with acceptable goodness of fit. A 
correlated 2-state fit to t-range 2-10 was found to give the most reliable determination of the 
ground state mass. The statistical error was evaluated using bootstrap. A systematic error from 
the uncertainty in this fitting procedure was determined by considering a range of other fits (by 
varying the i-range, fixing the excited state mass excitation to 1.5 lattice units, using a 1-state 
fit only, using an uncorrelated fit, fitting only subsets of the data). In each case this systematic 
error was found to be not larger than the statistical error quoted. 

Our lattice results for the 1~ + , + ~ and 2 +_ exotic hybrid mesons are given in the Table. 
Note that correlation studied in the last row in the Table can have contributions from a non- 
exotic 63 meson as well as the exotic 2" 1 , so it is only included for completeness. Our results are 
consistent with the assumption that the 1 h state is the lightest. Taking account of statistical 
correlations, we find that the + ~ state is heavier than the 1~ + state with a significance of 1 
standard deviation. 

As a cross check we also studied these states using a loop of colour flux which starts and 
finishes at the same site to model the gluonic excitation at the source. This closed gluon 
loop with non-trivial colour has the advantage of creating hybrid states which acn be studied 
without any additional quark propagators beyond those usually used in hadron spectroscopy. 
The disadvantage is that it creates hybrid states with the quark and antiquark at the same 
location and this may be an untypical component of the wavefunction. For the J PC = 1~ + 
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Figure 1: The lattice effective mass for the J = 1 h hybrid meson versus time separation 
t. The source used was a U-shaped path of size 6x6, while the sinks were combinations of 
U-shaped paths of size 6 x 6 (o), 3 x 3 (□) and lxl (fancy cross). We also used closed colour 
loops at the source of size 6x6 (*), 3x3 (octagon) and 1 x 1 (x) with the above sinks. The 
simultaneous two state fit to all six correlations is shown. 
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state, we used a 'magnetic' colour loop (L = 1 H ) and, as shown in the Figure, the effective 
mass from measuring 294 configurations is consistent with that obtained from U-shaped sources. 
From a fit to the closed-loop source correlations only, we obtained Ma = 1.12(8). For the 
other exotic states, the loop was taken to be in a L PC = 1 ^ representation and the resulting 
correlations were very noisy compared to those obtained from U-shaped sources. 

We express our lattice results for the hybrid masses as a ratio to the vector meson mass we 
find on the same lattices, namely Mya = 0.54(1). Since for the hybrid mesons, we find Ma ~ 1, 
we must address possible order a errors which do not cancel in this mass ratio as a is decreased 
to the continuum limit. The SW-clover fermionic action has been used with the intention of 
reducing such errors. Indeed a study |l(| varying a, shows that dimensionless mass ratios are 
consistent with constant over a range of lattice spacings including that used here, unlike the 
case when the Wilson fermionic action is used. It is difficult to estimate the systematic error 
from this extrapolation to the continuum limit of our results without repeating the calculation 
at smaller a. It is plausible that it is no larger than the statistical errors we quote. 

Another potential source of error is that we use a finite spatial volume of extent approximately 
1.6 fm. Investigations using larger volumes have shown that this is large enough for conventional 
mesons. For glueballs, finite size effects are known to be unimportant for spatial extent greater 
than 1 fm. The static quark analysis of hybrid mesons, however, finds that the effective potential 
is rather flat versus inter-quark separation so that the wave function will be spatially extended || 
Q . This implies that hybrid mesons may be significantly more extended than conventional quark 
model mesons or glueballs. Without repeating the hybrid mass calculation for larger volumes, 
we are unable to quantify the magnitude of finite size effects. 

For computational reasons, we have used a quark mass which is heavier than the strange 
quark mass in order to improve the signal to noise ratio. Ideally one should study hybrid 
mesons at a range of quark masses and interpolate to the experimentally relevant cases. Here 
we use instead a model assumption to make the small correction to our quark mass to bring 
it to the strange quark mass. This is not without some uncertainty as the model is based 
on experimental data whereas the lattice results apply to the quenched approximation - see 
ref. ]ll| for a discussion of this problem. In the quenched approximation, there will be magic 
mixing (i.e. no mixing) and the ss pseudoscalar will be at 0.68 GeV: thus M^/M^tg) = 1.5. We 
find My/Mp = 1.31(1) which implies that our quark mass is somewhat heavier than the strange 
quark. To correct for this, the simplest assumption, which is string inspired, is that for two types 
of quark q\ and q"2 : M(qi) 2 — M{q2) 2 is the same for hybrid, vector and pseudoscalar mesons. This 
implies that if our lattice ratio is M/My = x then for s-quarks M s jM^ = ^Jx 2 + 0.33(x 2 — 1). 
This correction is shown in the Table. 

Using the experimental cj> meson mass to set the scale, our result for the lightest ss hybrid 
meson is 1.99(13) Gev. Taking account of some of the possible systematic errors described above, 
we round this prediction to 2.0(2) GeV. We can use the model described above to estimate that 
the corresponding non-strange hybrid meson would be 120 MeV lighter than the ss state. 
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It is interesting to compare our result with model estimates. There have been two classes of 
model which have been proposed to describe the hybrid spectrum. Constituent gluon models, of 
which the bag model |[^] is the prototype, have a lowest gluonic excitation which is 'magnetic' 
with L PC = This implies that the J PC = 1"+ spin exotic hybrid will be lightest in 

constituent models. The spectrum of these 'constituent gluons' has been studied on a lattice |l4| 
which confirms that the magnetic gluon is indeed the lightest such excitation. The other class 
of model treats the gluonic excitation as an excited state of a string or fluxtube, as motivated 
by the lattice description Q of hybrid mesons with heavy quarks. A continuum version of 
this fluxtube model |l|] has been used to give predictions for hybrid mesons containing light 
quarks. These levels are expected to be relatively heavy: about 2.1 GeV for the ss system. 
This energy value comes basically from the string model which gives excited levels with energies 
nir/R higher for integer n. Then R, which is the interquark separation, can be estimated from 
standard quark models. This fluxtube model leads to degenerate hybrid levels which include 
spin-exotic J PC = 1 '", 0" 1 , and 2 H . A reasonable compromise from these two models is 
that the fluxtube model would be a good approximation to the hybrid spectrum but that the 
degeneracy of the spin-exotic levels would be broken in the direction predicted by the constituent 
aproach. This is exactly what our results indicate. 

A study of light-quark hybrid mesons needs input from the theoretical mass values in the 
quenched approximation to guide experimental searches and to provide input to models for the 
mixing between the qq hybrid states and qqqq or two meson states. Here we provide the first 
results of such a study from first principles. We focus on the spin-exotic states to avoid mixing 
with qq mesons. For these, the J FC = 1 ^ state is consistent with being the lightest. The 
mass value that we obtain is 2.0(2) GeV for the ss state. It is plausible that this error estimate 
includes the systematic errors but these systematic errors can be fully controlled by extending 
the lattice study to larger lattices, smaller lattice spacings and several quark masses. Our result 
for the mass spectrum of spin-exotic hybrid mesons is comparable to that suggested by fluxtube 
models. 

We acknowledge support from EPSRC grant GR/K/41663 to the UKQCD Collaboration. 
We thank Stephen Pickles and David Richards for their help in setting up the inversion codes 
for the shifted local propagators. 
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